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^N I Abstract. We show that in two Higgs doublet models at tree-level the potential minimum 

preserving electric charge and CP symmetries, when it exists, is the global one. Furthermore, 
C^ ■ we derived a very simple condition, involving only the coefficients of the quartic terms of the 

^ ', potential, that guarantees spontaneous CP breaking. 

o 
o 



In the Standard Model (SM) of the electroweak interactions the existence of the scalar Higgs 
Q^l doublet is a fundamental piece of the theory. Through it SU{2)w x U{1)y gauge invariance is 

jis ■ broken, the W^ and Z^ bosons and the fermions acquire their masses and the renormalisability 



C^ 



D , of the theory is preserved. Despite its importance, the scalar sector of the SM has not yet been 

'^ ' directly tested and there is considerable interest in studying its extensions. The simplest of 

those extensions is the two Higgs doublet model (2HDM). One of the main reasons of interest 
k> , in this class of models is the possibility of having spontaneous CP violation ^, thus helping 

J^ I to solve the baryogenesis problem P] (for a review, see 0i). One problem of these models, 

though, is their immense parameter space - the most general potential that preserves the SM 
gauge group and does not explicitly break CP has ten independent parameters, and little is 
known of their allowed values. A similar difficulty affiicts Supersymmetric models, where the 
parameter space is generally larger. One idea that has been applied to Supersymmetric theories 
to restrict their allowed parameter space is to use charge and colour breaking (CCB) bounds. If 
a given combination of parameters causes the appearance in the potential of a minimum where 
charged/coloured fields have vacuum expectation values (vevs), then that combination should 
be rejected. This appealing idea was introduced by Frere et al 4- and applied, in numerous 
papers, to several supersymmetric theories ^. Phenomenological analysis of supersymmetric 
Higgs masses use this tool to increase the models' predictive power 0|. It is therefore of interest 
to apply similar techniques to the 2HDM and try to limit its parameter space. The scalars of 
this theory have no colour quantum numbers but there are charged fields so charge breaking 
(CB) extrema are in principle possible. Recent work in 2HDM ^^, for instance, assumed that 
the choice of parameters made was such that the scalar potential respected the U{l)em gauge 
symmetry. 
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In reference jH] it was shown that to make sure that there were no stationary points corre- 
sponding to charge or CP spontaneous breaking, one had to restrict the parameter space to eight 
independent parameters. This leads to two independent potentials, stable under renormalisation 
because they are protected by a Z2 or U{1) global symmetries. It is interesting to stress that 
these are the usual symmetries introduced to prevent flavour changing neutral currents. In this 
letter we go back to the study of the tree-level vacuum of the most general 2HDM potential 
without explicit CP violation, with ten parameters. We prove that for this potential the min- 
imum that preserves CP and C/(l)em (from now on named the "Normal minimum"), when it 
exists, is global. This is a very powerful result, in that it assures that the Normal vacuum is 
stable, it cannot tunnel to other minima. This letter is structured as follows: we will review 
the model and the Normal minimum in section ^ then proceed to look at the possibility of 
charge breaking stationary points, proving there can be no CB minima if the potential has a 
Normal minimum. Likewise, in section |21 we will prove the equivalent result for spontaneous 
CP breaking stationary points. In section 0] we will look at the mass matrices of the three types 
of stationary points. 



1 The Normal minimum 

We will be working with the most general two Higgs doublet model invariant under the gauge 
group SU{2)w X U{1)y that does not explicitly break CP, following the conventions of refer- 
ence 0. This model is built with two scalar doublets of hypercharge Y = 1, 
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9^1 + ^'/^2 

995 + iipj 



<l>2 



ip3 + i(p4 

<^6 + i^s 



(1) 



This numbering of the ip fields might seem odd, but it is the most convenient for the latter calcu- 
lation of the mass matrices. As was shown in jHI this potential has ten independent parameters 
and may be written as 



V 



=ai xi + 02 2:2 + asx^ + 
611 xl + 622 xl + 633 X3 + 644 X4 + 612 a;iX2 + bi3 X1X3 + 623 3;2a;3 
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with 
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X4 



Xl = \<^i\^ 

X2 = 1^2!^ 



(piifS + (p2fA + f5'P6 + y'TV'S 
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The Hj parameters have dimensions of squared mass, the bij parameters are dimensionless ^, the 
fields ipi are real functions. Let us introduce a new notation that will be extremely useful: we 
define a vector A and a square symmetric matrix B as 
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^ These parameters have well-known relations to the masses and couplings of the physical Higgs particles of 
the model, see, for instance, |7|. 



Defining the vector X = {xi , X2 , X3 , X4), we can rewrite the potential Q in the more compact 
form 



V 



\T 
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A^ X + -X' BX 



(5) 



The Normal minimum corresponds to (p^ = vi-, '■P& 
This gives, from the above definitions, xi 
the relevant minimisation conditions as 



vl, X2 



V2 and all the remainder ipi equal to zero. 
vit;2 and X4 = 0. We can write 



V^, X3 



dV 
dvi 
dV 
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dV dxi dV 8x3 

^ — — + 
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dV dx2 

0^ - — -r^ + 



3x3 dvi 
dV dxs 
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(6) 



Let us define the vector V' , with components V- = dV/dxi, evaluated at the minimum. From 
the above it is plain that 



v 



'V{ 



2i;ii;2 



-2t;it;2 




(7) 



Looking at the expressions for the two first components of V' it is obvious that, regardless of the 
values of f 1, V2 and Kj', VI and V2 have the same sign. From this point forward we will use X^ to 
designate the vector X evaluated at the minimum, that is, with components {v\ , f | , "Wi f 2 1 0). 
In this notation it is trivial to realize that Xjj V 
shows that we can write V' in matrix form as 



0. Direct analysis of the potential © also 



V' = A + BXn 



(8) 



The potential (jSJ is a sum of quadratic and quartic polynomials, let us call them p2 and P4; by 
performing the sum ^^ Vi (dV/dvi) it is very simple to show that the minimisation conditions 
imply 2p2 + ^Pi = at the minimum. As such the value of the potential at this stationary 
point - which we designate by Vn - may be written as: 



Vn 



1 



A' Xn 



Xjsf B Xm 



(9) 



We have been speaking of the Normal minimum but the conditions © only assure us that the 
potential has a stationary point. To ensure we are at a minimum we must analyse the second 
derivatives oiV - that is, the matrix of the squared scalar masses - and reject all combinations 
of parameters {aj , hjk} for which any of the non-zero eigenvalues are negative (this matrix has 
three zero eigenvalues corresponding to the Goldstone bosons, see section^. In particular we 
observe that the squared mass of the charged Higgs is given by M^± = V^ + ¥2- So the Normal 
minimum exists if y/ + V^' > 0. Since we have already shown that Vl and V2 have the same 
sign this tells us that both quantities are positive. Therefore we obtain 



Vi 



and so we conclude that the quantity 
mental later on. 



Vj 

2i;iW2 



vl 



> 



(10) 



-K3'/(2t'it>2) is positive. This conclusion will be funda- 



2 Charge breaking 



The potential Q has only three types of non-trivial stationary points 001: the Normal one; a 
charge-breaking stationary point, where one of the charged fields (p has a non-zero vev; and a 
CP-breaking stationary point where one of the imaginary neutral component fields has a non- 
zero vev. In this section we will deal with charge breaking (CB), specifically, the configuration 
where the fields that have vevs are 935 = ii^, (/^g = v'2 and ip^ = a. The last vev breaks 
charge conservation and would give mass to the photon. We can calculate the derivatives of the 
potential with respect to f^, V2 and a explicitly and arrive at the conclusion that there is always 
a solution for a = 0, which corresponds to the Normal extremum. It is however easier to deal 
with the results if one uses the matrix notation. We define the vector Y to have components 
(xi , X2 , a^s , X4) evaluated at the CB stationary point, that is. 



Y 



,,'2 
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Av'2 




(11) 



Unlike the Normal case we now have three independent variables, this allows us to write the 
stationarity conditions as 



dV 
dX 



<^ A + BY 







X=Y 



Hence, as long as the matrix B is invertible, the solution Y is such that 

Y = -B'^A . 



(12) 



(13) 



We observe that the CB vevs are given by a linear equation, which means that this solution, 
if it exists, is unique. This is unlike the Normal case, where the minimisation conditions © 
lead to a system of two coupled cubic equations and as such can in principle produce multiple 
solutions. Notice however that the CB solution does not always exist even if B is invertible - the 
first two components of the vector Y must necessarily be positive, and not all choices of A and 
B matrices will give such a result in eq. H13|). Finally, the same reasoning that led us to eq. Q 
can be applied to this stationary point and we may write the value of the potential, Vqb, as 



VcB 



1 



a'^y 



1 r 
-Y^ BY . 



(14) 



We now have all the ingredients necessary to show that, if the Normal minimum exists, it is 
always deeper than the CB stationary point. We look again at equation Q and use equation H12|) 
to write 

V' = -BY + BXn . (15) 



Remembering that XTr V 



and thus 



we obtain 



Xj^BY + xIbXn 



Xjv BY = X^ B X]\i = — 2 Vn , 
the last step arising from equation Q. We now calculate the quantity y^y. 



(16) 

(17) 



Y^V 



Y^ BY + Y'^ BXn 



(18) 



Because the matrix B is symmetric we have Y B X^ = Xj^ B Y 
tion (UH), Y'^ BY = -2Vcb, so 



VcB — Vn 



■Y^V 



2 V/v and, from equa- 
(19) 



The product y^ V can be exphcitly written as (from ((Tj) and Hill) ) 
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At the end of section ^ we have demonstrated that the quantity — VJj'/(2i)it;2) is positive so we 
conclude that 

VcB - Vn > ^ , (21) 

which is to say, the Normal minimum is always deeper than the CB stationary point. So we 
conclude that the Normal minimum is perfectly stable and can never tunnel to a charge-breaking 
vacuum. In fact, as we will see in sectional the matrix B determines the character of the CB 
stationary point - if it is positive definite, so is the matrix of the squared masses at the CB 
stationary point. Multiplying both sides of eq. H15|) by B~^ first and then by V , it is simple 
to rewrite Y V' as 

Y'^V' = -V''^B'^V' . (22) 

We have shown that when the Normal minimum exists V' B^^ V' is negative. This implies 
that the matrix B^^ is not positive definite. So B is also not positive definite. We will prove 
in the appendix that its first entry, 2bii, is necessarily positive to prevent the potential from 
being unbounded from below. Hence B cannot be negative definite either. Therefore, the CB 
stationary point is necessarily a saddle point. 



3 CP breaking 



Besides the Normal minimum or the CB stationary point, we have another possible stationary 
point, one that spontaneously breaks CP conservation. In this case the fields which have non- 
zero vevs are c/35 = v'l, ip^ = v'2 and ipj = 5 - this last one breaks CP. The variables x,:, at 



this stationary point, are xi 



,,'2 



+ 5\ X2 



/2 
V2 , X3 



f]^f2 and X4 



f 2 6- We see that 



xi X2 — xi and as such the potential for this field configuration can be written as 



V =ai xi + 02 X2 + 03X3 + 

„2 , r,__2 , /.__ L \^2 



felia^l + 622^:2 + (^33 - b44)x^ + (612 + 644)2;iX2 + bisXiX^ + 623X2X3 

Defining the vector Z and the square, symmetric matrix Bcp to be 
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(24) 



we see that the potential at this stationary point may be written in terms of {v[ , V2 , 6} as 

1 



VcP = A^Z + -Z'^BcpZ . 



(25) 



Then we are exactly in the conditions of the previous section: the CP breaking vevs are given 
by the equation 



Z 



BcP^ 



the value of the potential at this stationary point is given by 



V, 
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— Z Bnp Z 
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and repeating the steps of the previous section we find 

Vcp — Vn = —- 
Calculating Z^ V explicitly we find 



'-rV'^B^lV' 



-z^V 

2 



Z^V 
and so we conclude that 
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(26) 



(27) 



(28) 



(29) 



(30) 



Vcp — Vn 

That is, once we have found a Normal minimum, it is always deeper than the CP-breaking 
stationary point. Again, we have proved that Bqp is not positive definite and, since its first 
entry is positive (see the appendix), it cannot be negative definite. As before Bcp determines 
the nature of the squared mass matrix at this stationary point (see section |1}, which means that 
if a Normal minimum exists, the CP stationary point is necessarily a saddle point. 



4 Mass matrices 



To determine the nature of the stationary points one must analyse the second derivatives of the 
potential, which is to say, the scalar squared mass matrices. They are given by 

d'^V dV d'^xi d'^V dxi dxr 



difidipj 
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dxi difidipj dxidxm dipi dipj 



(31) 



Extending the notation of section ^ we define V^' = dV/dxi evaluated at any of the three 
stationary points that we are considering. The first term in the above equation is written as an 
8x8 matrix of the form 



[Ml] = 
where M^^ and M^2 ^^^ 4x4 matrices given by 
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(32) 



(33) 



In the second term of eq. (|3H) the derivative d'^V/dxidxm is clearly the matrix element Bim of 



the matrix B defined in eq. 
we call C, given by 



[C] 



\. As for the derivatives dxi/d(fj, they form a 4 x 8 matrix which 
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This matrix is evaluated at each of the different stationary points which is why only the fields 
{ip3 , c/95 , (/9e , (^7} appear, the rest are always zero at the stationary points. Then the scalar 
squared mass matrix becomes 



[m2] = - {[Mf] + C^BC) 



(35) 



where the factor of 1/2 in the left hand-side is due to the fields ^p being real. Let us now analyse 
this matrix for each of the three cases. 

• The Normal stationary point 

Prom section 121 we see that Ki = and V^ = AV^ V2. In this case only ip^ = vi and ipQ = V2 
are non zero, so that the first four columns of the matrix C are zeros. As a consequence the 
matrix C B C is block diagonal, with the first diagonal 4x4 block composed of zeros, that is, 



C^BC 



0' 
B' 



(36) 



where 5' is a 4 x 4 matrix. In our notation this means that the matrix Ml^/2 is now the mass 
matrix of the charged Higgs, and it is very simple to see that it has two zero eigenvalues and 
a doubly degenerate eigenvalue given by ¥{ + ¥2- The matrices M^2 and B' are also block- 
diagonal, with two 2x2 non-zero matrices along the diagonal. One of these is the mass matrix 
of the pseudo-scalar sector, given by 



Vl + 644 V2 V3 — 644 vi f 2 

Vg — 644 Vl f 2 V2 + 644 V2 



(37) 



This matrix has a zero eigenvalue and another one equal to V^ + V2 + ^44 (ff + v^)- The mass 
matrix for the CP-even scalar sector is given by 



WIh] 



Vi + Hi Vl + Hz 
Vi + Hz Vi + H2 



(38) 



with 



Hi=A hii vf + 2 613 Vl V2 + &33 V2 

H2 = 4: 622 «2 + 2 &23 Vl V2 + 633 vl 

H3 = (2612 + 633)^1^^2 + ^13^1 + b23V2 



(39) 
(40) 



• The CB stationary point 

The mass matrix in this case is very simple given that every V^' is zero, which means M^ = 0. 
The matrix C now has only one column of zeros but it is very easy to see that three other columns 
are linearly dependent. With judicious operations on the lines and columns of the matrix C we 
manage to set to zero its first four columns and therefore write 
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ro 







/T 



C" B C 



(41) 



where C" is a 4 x 4 matrix. It is then obvious that M^ has four zero eigenvalues, exactly those 
we would expect having broken the U{l)em symmetry. The non-zero eigenvalues are then those 
of the matrix 



[M< 



CB\ 



- C'^ B C 
2 



(42) 



and this expression implies that [M^^] is positive definite if and only if B is positive definite. 
If B is negative or semi-positive definite [-/W^^fi] ^^ likewise defined. This justifies our assertion 
in section 121 that the CB stationary points are saddle points. 

With further operations on lines and columns it is still possible to simplify [M^^]: one of its 
eigenvalues is 644 {v'l + f 2 + a^), the remaining three are those of the following 3x3 matrix, 



with 
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CB33 = bi3v'iv'2 + Uiiv'i^ . 
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The CP stationary point 



Again the matrix C has the first four columns equal to zero, which means that the matrix B 
only has an impact on the lower right 4x4 corner of [M^] . From the expressions of section |31 it 
is easy to obtain V{ = — 644 v'2 , V2 = — 644 {v'l + 6"^), V^ = 2 644 v'l v'2 and V'l = —2 644 v'2 5. 
Once again the matrix Mf is the squared mass matrix of the charged Higgs sector, it has two 
zero eigenvalues and a doubly degenerate non-zero one, given by —644(1;']^ + ^2 + f^^)- The 
6 vev now causes mixing between the CP-even and odd scalar fields so that we are left with a 
4x4 symmetric matrix, M^p, whose entries are 

Ml,pil,l)=4biiv'i^ + 2bi3v'iv'2 + (633 - M4' 

Mcpil, 2) = 613 v'l + (633 -h 644 -h 2 612) v'l v'2 + 623 v'2 

Ml:p{l,3) = {bisv'2 + Uiiv'i)6 

M2p(l,4) = [bisv'i + (633 - b^4)v'2]5 

M^p(2,2) = (633 - bM)v'i^ + 2 623 4 ^^2 + 4 622 4' 

M^p(2,3) = [bi3v'i + 2(612 + ^44)^^2]'^ 

M2p(2,4) = [623^^2 + fe - b^4)v'i]6 



M^p(3,3)=4 6n52 

M^P (3, 4) =6i3 (^2 

M2p(4,4) = (633 - bu)S^ . (45) 

This matrix has one zero eigenvalue which gives a total of three Goldstone bosons, as was to be 
expected. The most interesting aspect of the CP mass matrix lies in the value of the squared 
charged mass, here proportional to —644, whereas in the CB case it was proportional to +644. 
This means that we can never have, for the same choice of 644, a CB minimum and a CP one, 
as was concluded in ref. [S] • If we define the matrix C" to be the restriction of matrix C to its 
last four columns it is trivial to show that 

[M^p] = \c"^BcpC" . (46) 

So, just like the CB case, the Bcp matrix determines the nature of the stationary point: if it is 
positive definite - and 644 < 0, due to the charged Higgs eigenvalue - the stationary point is a 
minimum. If not, it is necessarily a saddle point. 



5 Conclusions 

In this paper we have shown a remarkable result: that the two Higgs doublet model is "protected" 
against electric charge or CP spontaneous breaking. In other words, if the model has a minimum 
preserving C/(l)em and CP, that minimum is global. In this way, there is absolutely no possibility 
of tunneling to deeper minima, and, for instance, the masslessness of the photon is guaranteed 
in these models. We also obtained a simple relation telling us that the only case where charge 
breaking might occur is when the matrix B is positive definite. This situation implies that no 
Normal minima exists. For the future, when one is scanning the parameter space of the 2HDM 
one can safely exclude, from the beginning, the combinations of bij parameters that give rise to 
a positive defined B. 

Charge breaking would be disastrous but there is considerable interest, from cosmologists 
to particle physicists, in models with the possibility of spontaneous CP violation. We have 
determined that this cannot happen for those ranges of parameters that lead to Normal minima. 
However, we have also established a very precise condition for spontaneous CP breaking to occur: 
CP is spontaneously broken if and only V' B^p V > Q^. In these circumstances the 2HDM no 
longer has a Normal minimum. A simpler condition for spontaneous CP violation is to demand 
that the matrix Bcp be positive definite - this condition, together with 644 < 0, guarantees 
that the CP stationary point, when it exists, is a minimum. We can also guarantee that if we 
find a CP minimum, that too is safe against charge breaking. This arises from the fact that a 
CP minimum requires 644 < which automatically implies the matrix B is not positive defined. 
Then we will have Vcb — Vn > and Vcp — Vn < 0. Therefore, Vcb — Vcp > and the 
CP minimum is safe against charge breaking. 

Let us also stress that our conclusions are absolutely general, independent of particular 
values of the parameters of the theory, obviously. They hold for any of the more restricted 
models considered in ref. 8 . It is simple to recover the conditions presented in that reference 
to avoid CP minima by analysing the matrix Bqp- We remark that the Higgs potential of the 
Supersymmetric Standard Model (SSM) is also included in the potentials we studied - in fact, 
it corresponds to the case 611 = 622 = — &12/2 = Af^/(2t;^), 633 = 644 = IMy^/v"^ and the 

^In reference ^ CP violating quantities involving only the Higgs sector were derived in models with explicit 
CP violation. 



remaining b parameters set to zero, following the conventions of ref. jEIl- So we could conclude 
that at tree-level, the Supersymmetric Higgs potential is safe against charge of CP violation, 
though this would not preclude charge, colour or CP breaking arising from other scalar fields 
present in those models. However, we must be cautious: it has been shown ^J that one-loop 
contributions to the minimisation of the potential have an enormous impact on charge breaking 
bounds in Supersymmetric models. Also, it was recently shown pl2 that unless one performs 
a full one-loop calculation (both for the potential and the vevs, in both the CB potential and 
the "normal" one) the bounds one obtains can be overestimated. Therefore, we urge caution in 
applying these conclusions to the SSM. Nevertheless one would expect the one-loop contributions 
to be much less important in the non-supersymmetric 2HDM due to the much smaller particle 
content of the latter theory. 
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Appendix: bounds on the 6-paranieters 

Requiring that the potential ^ be bounded from below in all possible directions imposes some 
interesting bounds on the b parameters. These bounds are obtained studying particular field 
directions and requiring that the quartic terms' limit as the fields go to infinity be positive or at 
most zero. For instance, we can choose a direction such that all variables x except xi are zero, 
and xi ^ oo (for instance, by choosing 991 — > oo and all remaining ipi equal to zero). Then along 
this direction the potential goes to infinity as buipf, and if we want this limit to be positive 
or zero, we must demand that bn > 0. Likewise, by choosing ip^ ^>- 00 and all others zero, for 
instance, we would obtain the condition 622 > 0. We cannot obtain a similar condition for 633 
and 644 because we can never find a field direction with ^3^4 -^ 00 without having xi^2 diverging 
as well. Let us now consider the direction (/?i, (^5 — > 00 and all others set to zero. Only xi and X2 
are non-zero and the potential is reduced to the terms buipf + 622 ^5 + &12 ^1 ^\- By choosing 
polar coordinates such that (p\ = r cos 9 and ip2 = r sinO, the potential will always be greater 
or equal to zero at infinity along this direction if, for any < 9 < '7r/2, 

611 cos^ e + 622 sin^ 9 + 612 sin 6* cos 6* > . (47) 

A simple minimisation in 9 shows that this occurs as long as 



612 > - 2^611 622 • (48) 

Likewise, if we choose the direction (/?!, c/?4 — > 00 and again use polar coordinates, the bounded- 
ness from below condition translates into 

611 cos2 6' + 622 sin^^* + (612 + 644)sin6'cos6l > (49) 

and so the condition we obtain is similar to the previous one, 

bi2 + 644 > -2V611622 . (50) 

With the direction 931 , (/93 — > 00 we can study what happens with the b^i parameters, the condi- 
tion we obtain is (making ipi = r cos 9 and (^3 = r sin 9 so that now there are no restrictions 
on the value of 9) 

bucos'^9 + b22sm^9 + (612 + 633) sin^ 6'cos^ 6* + 613 cos^ 6* sin + 623 sin^ 6^ cos 6* > . (51) 



Since this condition has to hold for any value of 9 we can make the change 6^—0 and obtain 
|6i3 cos^6'sin(9 + 623 sin3 6'cos(9| < 6iicos^(9 + b22sm'^0 + (612 + 633) sin^ 6'cos^ 6* . (52) 
So, we conclude that, for any 6, we have 

(611 cos'' 6* + 622 sin^^* + (612 + 633)sin2 6'cos2 6l > (53) 

and we obtain a condition similar to 1)481) and (|5U() . 



bi2 + bss > -2y/bubT2 ■ (54) 

Finally, making 6 = tt/4 in eq. H52() . we obtain a more manageable bound on 613 and 623, 

1^13 + &23I < ^11 + ^22 + bi2 + ^33 • (55) 
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